Isotope effect and the role of phonon in the Fe-based superconductors 
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We studied the isotope effect of phonon in the Fe-based superconductors using a phenomenological two band 
model for the sign-changing s-wave (±s-wave) state. Within this mean-field model, we showed that the large 
isotope effect is not inconsistent with the ±s-wave pairing state and its high transition temperature. In principle, 
a large phonon isotope coefficient a implies a large phonon coupling constant. However, the asymmetric density 
of states (DOS) between two bands substantially enhances the value of a, so that a moderate value of the phonon 
coupling constant (X p j, ~ 0.4) can produce a very large value of a (~ 0.4) as well as a high transition temperature 
together with an antiferromagnet (AFM)-induced interaction. 

PACS numbers: 74.20,74.20-z,74.50 



Recent discovery of the Fe-based superconductors by 
Kamihara et al. [01l2l]> nas greatly spurred the research activ- 
ity of unconventional superconductivity. Regarding the pair- 
ing mechanism and symmetry of theses new superconducting 
(SC) materials, there are already numerous experimental and 
theoretical investigations. Since the first theoretical proposal 
of the ±s-wave state by Mazin et al., [3] as a best pairing 
state in the Fe-based superconductors (SCs), several subse- 
quent theoretical studies ansa supported this idea. Ex- 
periments such as ARPES [8] and penetration depth measure- 
ments [9] unanimously indicate a full gap around the Fermi 
surfaces (FSs) consistent with a s-wave gap state. Then the 
nuclear spin-lattice relaxation rate 1 /T\ 1 10], which is seem- 
ingly consistent with a nodal gap state such as a d-wave gap, 
provided a strong evidence of the sign-changing nature of the 
gaps on different bands, so that actually strengthened the case 
of the ±s-wave state lH [TTl [Hi . 



As to the pairing glue, most researchers at the moment tend 
to believe an electronic origin rather than a phonon origin par- 
ing flBEBGH. In particular, an antiferromagnetic (AFM) 
correlation induced interaction appears to be the most natural 
pairing glue in view of the common SDW instability at around 
^ 150 K and the overall phase diagram with doping in this se- 
ries of Fe pnictides [ 13J . It is shown by several authors that 
the AFM-induced potential, when combined with the unique 
band structure (or the FS topology) of the Fe pnictides, nat- 
urally leads to the ±s-wave state as the best pairing solution 



However, Liu et al. [ 14] recently measured isotope effect on 
T c with a substitution of 56 Fe by 54 Fe in Bai_ x K x Fe2As2, and 
reported a unexpectedly large isotope coefficient a (~ 0.4). 
This observation, if confirmed, is drastically perpendicular to 
the current line of thought for the pairing mechanism of the 
Fe-based SCs. Theoretical investigations about the electron- 
phonon coupling in these materials are yet only a few. Boeri et 
al. lfl5ll calculated a very weak electron-phonon coupling con- 
stant (X p h < 0.2) - the value averaged over the Brillouin Zone 
(BZ). Eschrig I1 1611 argued that the proper electron-phonon 
coupling is not the averaged coupling constant but the one of 
a particular phonon mode, i.e., in-plane Fe-breathing mode, 



which may have a large electron-phonon coupling constant. 

Besides theoretical investigations, the experimental fact is 
that all Fe pnictides, either /?FeAsO (R= La, Ce, and Nd) or 
AKFe2As2 (A=Ba and Sr) compounds, display the structural 
instability from tetragonal to orthorhombic symmetry and it 
always occurs at temperatures very near the SDW transition 
temperatures 01711 . It implies, at least, two things: (1) the lat- 
tice degrees of freedom (structural instability) and the spin de- 
grees of freedom (magnetic instability) are not independent; 
(2) this structural instability is not a Jahn-Teller-type instabil- 
ity of the Fe 3d-electrons but is closely related with the metal- 
licity of the Fe 3d-electrons 01 811 . Therefore, the issue of the 
electron-phonon coupling in the Fe pnictides needs to be fur- 
ther investigated. 

With this motivation, in this paper, we studied the isotope 
effect and the role of phonon using a phenomenological model 
for the ±s-wave pairing state. The details of the model can be 
found in Ref.|5] and here we briefly sketch the essential in- 
gredients for our purpose. The model consists of two bands: 
one hole band centered around Y point and one electron band 
centered around M point in the reduced Brillouin Zone (BZ) 
scheme, and has a phenomenological pairing interaction, in- 
duced from an AFM correlation, hence peaking around (tc,tc) 
momentum exchange. In this paper, we add a phonon inter- 
action to this model. Because our purpose of this paper is 
to study the isotope effect of phonon when the total interac- 
tion gives rise to the ±s-wave pairing state, we assume only 
a general condition of the phonon interaction and vary the ba- 
sic parameters of the phonon interaction such as the coupling 
strength A. p h and the characteristic phonon frequency co p h. 

The Hamiltonian is written as 
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where Z ne {k) are the dispersions of the hole band and elec- 
tron bands, respectively, representing two main bands in the 
Fe pnictides. The details of the dispersions are not important 
for our purpose except the density of states (DOS) of each 
band, N n (hole band) and N e (electron band), respectively. h' ka 
and e T ka are the electron creation operators on the hole and 
the electron bands, respectively. As mentioned previously, 
Vafm(&,£ ) is the AFM-induced pairing potential, which is 
all repulsive in momentum space, and V v h(k,k ) is the phonon 
interaction, which is all attractive in momentum space. 

The minimum characteristics of the interactions to promote 
the ±s-wave gap solution are: Vafm(^,£ ), peaking around 
(7t,7c) momentum exchange, should have a stronger inter- 
band interaction than the intraband one; and Vph(k,k ), be- 
ing stronger for small momentum exchange, should have a 
stronger intraband interaction than the interband one. The lat- 
ter requirement for the phonon interaction is already included 
in the Hamiltonian by not including the interband terms like 
Vph(k,k )*j[t*!jfcjA it /jA_ jfc /.|, and its hermitian conjugate. This 
assumption of the phonon interaction is indeed the property 
of the main phonons in Fe pnictides II 1 511 . If the phonon inter- 
action were absolutely momentum independent, it would have 
null effect for the ±s-wave pairing. 

For simplicity of the analysis but without loss of 
generality, we only need the FS-averaged interactions: 
for the AFM-induced interactions such as V^fm =<< 

VAFM(kh,k e ) »k k ,k e > VaFM =« VaFm(*A,*a) »k,„k' r etc - 



and similarly for the phonon interactions such as v££ 

V ph {k h ,k' h ) » kh ,kf h andy ph = - «Vph(k e ,kf e ) »k e ,k>- No- 
tice that, in these definitions, we absorbed the signs of the 
interactions and therefore all Vj^ M and are positive val- 
ues. Also assuming the ±s-wave solution we fix the signs of 
the s-wave gaps as A/, = \A„ | on the hole band and A e = — \A e \ 
on the electron band, respectively. The coupled r c -equations 
are written as 



= - << 



A/, 



A, = 



{V& M N eX AFM ]A e , 
[V& M N eX AFM -V^N eX ? h ]A e 

-[V& M N hX AFM ]A h . 



(2) 



(3) 



The pair susceptibilities X AFM and % ph at T = T c are defined 



as 
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where COafm an d <*>ph are the cutoff frequencies of the AFM 
fluctuations and phonon, respectively. The second expres- 
sion is the well-known BCS approximation valid only when 
<X>AFM,ph 3> T c , otherwise the first expression should be nu- 
merically calculated. Equations (2) and (3) with Eq.(4) con- 
stitute the r f -equation of the ± s-wave state of the two band 
model, including phonon pairing interaction as well as the 
AFM-induced pairing interaction. Before we show the nu- 
merical results we can analyze a simple case and gain general 
insight about the model. 

In the case (we call it "symmetric case") that N n = N e , and 
V^fm = ^afm and Vg = V** (V^ M = V^ M is always true), 
we define the dimensionless pairing constants as follows, 
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Then, the r t -equations are simplified as 



(5) 
(6) 
(7) 



A*(l+3^X A ™-VX ph ) = fflx^, (8) 
A,(l+«X AFM -W h ) = ^fmX AFM A„ 

The above T c -equations can be solved analytically with the 
BCS approximation of the pair susceptibilities X AFM and % ph , 
the second expression of Eq.(4), as 



^■ = 1.14cC M -c<-exp 



where 
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(10) 

(11) 
(12) 



Equations (9) and (10) show that the AFM interaction and 
the phonon interaction is additive in the exponential form so 
that T c can be dramatically boosted even with a small value of 
X p h I19I1 . We can also easily read the phonon isotope coeffi- 
cient a from Eq.(9) as 



a : 



1 d\nT c 

2 aflncOph 
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(13) 



This result conforms with a physical insight, i.e., a large a 
value means a relatively large phonon coupling compared to 
the total AFM pairing interaction X e // M "" = (?^fm - Kfm)- 
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FIG. 1: (Color online) Symmetric case (N^/N? = 1.0). (a) 
Isotope coefficients a as a function of the phonon cutoff fre- 
quency CDph/cOAFM for different phonon coupling constants Xph = 
0.1,0.2,0.3,0.4 and 0.8. (b) Normalized T c vs C0 ph /C0AFM for the 
same parameters. Symbols are numerical calculations and lines are 
the results of Eqs. (9) and (13). 
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FIG. 2: (Color online) Non-symmetric case (Nh/N e = 3.0). (a) 
Isotope coefficients a as a function of the phonon cutoff fre- 
quency cOph/cOAFM for different phonon coupling constants Xph = 
0.1,0.2,0.3 and 0.4. (b) Normalized T c vs co p h/cOAFM f° r the same 
parameters. Symbols are numerical calculations and lines are the 
results of Eqs. (9) and (13). 



However, Eq.(9) also shows that a large a value does not nec- 
essarily mean the large phonon contribution to the pairing en- 
ergetics when C0ph is much smaller than COafm- 

For the general cases where N/, ^ N e , and V^fm ^ ^afm an( l 
Vph ^ph (" non -sy mme tric case"), it is not possible to find an 
analytic solution of the r c -equation. However, an inspection 
suggests us to generalize Eq.(9) to the non-symmetric case. 
Equations (9) and (13) can be used as good approximations 
for T c and a with the definitions of the effective dimensionless 
coupling constants as follows, 

^AFM = d4) 



V = yJWeV$V$. (16) 

In the following, we will show numerical results directly ob- 
tained from Eqs. (2) and (3) and compare them with the an- 
alytic formulas (9) and (13) both for the symmetric and non- 
symmetric cases. 

Fig. 1(a) shows the isotope coefficient a of the symmetric 
case (Nh = N e ). Symbols are the numerical results of Eqs. (2) 
and (3), and lines are the results of Eq.(13). All the energy 
scales are normalized by the cutoff frequency of the AFM- 
fluctuations COafm an d we choose moderate strength of di- 
mensionless coupling constants of the AFM-induced potential 
as AS = N h V h ^ M = NeV&u = 0.6 and X^Sl = ^afm = 
iVeVfpM = 0.4, which yields the effective total magnetic in- 
teraction X e ^p M tot = 0.2. Fig. 1(a) shows that the numerically 
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calculated a (symbols) are in excellent agreement with the an- 
alytic results of Eq.(13) when cOph 3> T c . When this condition 
is not fulfilled, deviations occurs at low co p h where the BCS 
approximation of the pair susceptibility becomes poor. 

The results of Fig. 1(a) shows that a large phonon isotope 
coefficient arises when the phonon coupling strength is much 
stronger than the magnetic coupling strength and it is in ac- 
cord with a standard expectation. For example, in the sym- 
metric case, in order to obtain a w 0.4 (the reported value by 
Liu et al. ffl ). we need to have A p h —4x A^^/"' which cor- 
responds to the result of Aph = 0.8 (purple diamond symbols). 
If this is the case of reality, then the superconductivity of the 
Fe pnictides is a phonon driven SC and the AFM-fluctuations 
merely act as a tipping agent to introduce the % phase between 
two bands widely separated in the BZ. It is not an impossible 
scenario, but we first need to find an evidence of such a strong 
phonon coupling in Fe pnictides. If confirmed, the current 
viewpoint about the pairing mechanism of the Fe-based SCs 
should completely be changed. 

In Fig. 2 we show the results of the non-symmetric case 
(Nf, = 3N e ). We adjust the parameters to keep the dimen- 
sionless coupling constants A^^ = 0.6 and A^'J^ = 0.4, the 
same values as in Fig.l. Here A,^Jw" ra are calculated ac- 
cording to the generalized formulas Eqs. (14) and (15). The 
phonon coupling constants A p h = 0.1,0.2,0.3,0.4 are also the 
values according to Eq.(16). The behavior of a is very differ- 
ent compared to the symmetric case: (1) there are systematic 
deviations of the numerical (exact) results from the analytic 
(approximate) formula Eq.(13) and the numerically calculated 
values of a give substantially larger values than the ones of the 
analytic formula. As a result, a moderate value of A p h = 0.4 
can cause a large isotope coefficient a « 0.4; (2) with increas- 
ing (Mph, there is no saturation of a, in particular with large 
phonon couplings. In contrast, the T c values are rather similar 
to the symmetric case. 

This seemingly inconsistent behavior between T c and a 
can be most easily understood by comparison of the A p h = 
0.4 results between the symmetric case (Fig.l) and the non- 
symmetric case (Fig. 2). a in the symmetric case quickly sat- 
urates to its maximum value (red square symbols in Fig. 1(a)) 
but a in the non-symmetric case gradually increases with co p h 
( red square symbols in Fig. 2(a)). As a result, a large differ- 
ence of the a value between the two cases is possible when 
Wja/wafm — * 1 while maintaining only slightly enhanced T c 
(see red square symbols in Fig. 1(b) and in Fig. 2(b)). Slightly 
enhanced T c is due to the fact that the asymmetric phonon cou- 
plings on the hole band and electron band are slightly more 
advantageous for pairing than the equal phonon couplings on 

both bands when A ph = ^N h N e V$V£ is the same. 

In summary, using a minimal two band model with both the 
AFM-induced interaction and the phonon interaction, which 
together yields the ±s-wave gap pairing state, we studied the 
phonon isotope coefficient a for the symmetric case and the 
non-symmetric case. We confirmed that a large value of a 
indicates a large value of the phonon coupling constant A p h 



compared to the value of the AFM-induced coupling constant 
A-afm'"'- However, we found that the large asymmetric ra- 
tio (Nh/N e ) of DOSs between the hole and electron bands 
substantially enhances the a value compared to the symmet- 
ric band case. As a result, a relatively small value of the 
phonon coupling constant (say, A p h = 0.4 when X e /p M ""=0.2) 
can yield a very large isotope coefficient a=0.4. Regardless 
of the symmetric or non-symmetric cases, T c is strongly en- 
hanced by the addition of two pairing interactions A p h and 

^AFM m tne ex P onen tial form. A similar result was obtained 
for the high-r f cuprates SCs IU9I1 . The possibly important role 
of phonon(s) for the unconventional SCs discussed in this pa- 
per is a very plausible scenario, in general. Although it is 
studied based on a mean field theory, it is sufficient to demon- 
strate the principle. Further experimental study of the isotope 
effect in the Fe pnictide SCs is a pressing task. If the large 
isotope effect is confirmed, the importance of phonons for the 
unconventional SC pairing in the correlated materials should 
be renewed. 
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